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We suggest how to construct an effective low energy Hamiltonian via Monte Carlo starting from a given action. 
We test it by computing thermodynamical observables like average energy and specific heat for simple quantum 
systems. 
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1. INTRODUCTION 

Although most progress in lattice gauge the- 
ory has been achieved by use of the Lagrangian 
formulation, it is worthwhile to keep the Hamil- 
tonian formulation in mind. The potential useful- 
ness of the Hamiltonian lies in the following ar- 
eas: - Non-perturbative computation of scattering 
cross sections and decay amplitudes in hadronic 
systems. - Low-lying excited states of the hadronic 
spectrum and quantum chaos in such a system. 
- Hadron wave functions and hadron structure 
functions (for small xb and Q^). - Finite tem- 
perature and finite density in baryonic matter 
(quark-gluon plasma phase transition, neutron 
stars and cosmology). - Atomic physics: study 
of spectra and of quantum chaos. - Condensed 
matter, study of spin systems (computation of 
dynamical structure factors), and high Tc super- 
conductivity models (search for electron pair at- 
traction at very small energy). 

Here we suggest to construct an effective low 
energy Hamiltonian, denoted by Heff in the fol- 
lowing, via Monte Carlo. We suggest to apply the 
Monte Carlo method for the numerical computa- 
tion of matrix elements of the propagator, and 
second for a stochastic selection of states from a 
basis of Hilbert states. 



2. CONSTRUCTION OF H^ff 

Let us outline the method, by first assuming 
that we are given a finite set of Hilbert states, 
which is sufficiently large and suitably chosen 
such that low energy observables can be well re- 
produced. Later we will discuss the stochastic ba- 
sis selection. We start from a complete orthonor- 
mal basis of Hilbert states je.; >, i = 1, 2, 3 • • • and 
consider the matrix elements for a given fixed N. 
We consider a transition amplitudes in imaginary 
time 

M,,(r) e.le-^^/^le, >, i,j G 1, • • • , TV, (1) 

where H is the exact Hamiltonian. If H is Hermi- 
tian, then AI is unitarily equivalent to a diagonal 
matrix D such that 



M(T) = D{T) U. 



(2) 



By algebraic solution for D, we can construct ex- 
plicitly an effective Hamiltonian 



N 



k I 



(3) 



k=l 



where the \E'^^^ > are the eigenstates of D, and 
E'^^^ are the eigenvalues of log{D). One has the 
identities 

Ul^<e.\Elff>, D,{T) = e-^l"^'\ (4) 
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3. MATRIX ELEMENTS 

In order to be specific, we work in D = 1 di- 
mension. We choose basis states |ej > in posi- 
tion space by introducing a lattice with nodes Xi 
and define ei{x) (unnormahzed) by ei{x) = 1 if 
Xi < X < Xi+i, zero else. We choose Axi = 
Xi+i — Xi = const. The matrix elements read 

Mij{T) = I dy I dz 



J Xi J Xj 



'x] exp[—S[x\/K\ 



z.O 



(5) 



Here S denotes the Euclidean action for a given 
path C, 



S[C] = I dt ]-mx^ 
Jo 2 



V{x) 



(6) 



We suggest to compute the matrix elements 
Mij (T) from the action via Monte Carlo with im- 
portance sampling (Metropolis algorithm). In or- 
der to do so we write Mij{T) as a ratio of two 
integrals by splitting the action 



S = Sq + S\ 



dt —mx^ -\- 



dt (7) 



and expressing My by 



Mi^{T) = M\f{T) 



J Xi J Xj 



dz 



I 



[dx] exp[—Sv[x]/U] exp[—So[x]/h] 



z,0 



I / dy 



dz 



Ix] exp[— S'o[a;]/?i] 



z.O 



(8) 



Now we consider O = cxp[— 5*^/??,] as an observ- 
able and the matrix clement Mij is the expec- 
tation value of that observable. The matrix ele- 
ments M^j \ corresponding to the free action 5*0 
are known analytically (apart from an overall in- 
tegration), 



dy 



^3 + 1 



dz 



2wnT 



exp 



2nT 



(9) 



4. STOCHASTIC BASIS 



The success of lattice Euclidean lattice calcula- 
tions depends crucially on the efficiency of Monte 
Carlo with importance sampling when comput- 
ing the path integral. The path integral is rep- 
resented by a sum over "representative" configu- 
rations (corresponding to thcmodynamical equi- 
librium configurations from a Boltzmann distri- 
bution). Guided by that we suggest here to con- 
struct a basis of "representative" states such that 
the energies of the effective Hamiltonian follow 
(at least approximatively) a Boltzmann distribu- 
tion. The idea is to guide the basis selection 
by the quantum mechanical Euclidean transition 
amplitude. For the Free system it reads 



Geuci{x,T : y,o) = 



m 



2TTnT 

r m 
X exp— — — 



{x-yf]. (10) 



This function being positive for all x,y,T, can 
serve as a probability density. We choose y = 
and write 



P{x) = ^Geuci{x,T;0,0), 
Z = J dx Geuci{x,T;0,0). 



(11) 



Then we draw nodes Xi^ randomly, and construct 
basis states as box functions like those of the orig- 
inal nodes Xi . The idea is to work with a stochas- 
tic basis of dimension iVg// being much smaller 
than N, the dimension of the original basis. This 
should be feasible in particlular in many-body 
systems. In the free case P{x) is a Gaussian, 



Pix) 



1 , X 

exp 



27r(T 



2cr2J' 



(12) 



One can see that this implies a Gaussian distribu- 
tion for the wave numbers k and hence a Boltz- 
mann distribution for the energies E. In the pres- 
ence of a potential, a possible approximation is to 
write GeucI as a path integral retaining only the 
classical path. 
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Figure 1. Average energy of the harmonic osciUa- 
tor. Sohd hne and diamonds, respectively, repre- 
sent the exact analytical result, and that from the 
exact matrix elements for Ax — 1 and N = 20. 
Fig. taken from Ref.||l[]. 



Figure 2. 
oscillator 
Ref.§. 



Specific heat over ks of the harmonic 
Symbols as in Fig.[l]. Fig. taken from 



5. NUMERICAL RESULTS 

In order to test the Monte Carlo Hamiltonian 
we have considered the free system, the harmonic 
oscillator and other local potentials, where ei- 
ther analytic or precise approximative solutions 
are known. We have computed the energy spec- 
trum, wave functions and thermodynamic observ- 
ables like the partition function Z, the average en- 
ergy U and the specific heat C. Results for wave 
functions of the harmonic oscillator are shown in 
Ref. . Results for other potentials can be found 
in Refs. 1^,1). The thermodynamic observables 
are defined by 



uip) 

CiP) 



dlogZ 



= kB0 



dl3 ' 
logZ 



(13) 



where /? = {kBT)~^, T is the temperature, and 
we identify /3 with the imaginary time T by 
/3 — T/h. The partition function of the effective 
Hamiltonian is obtained from its spectrum. 



Zeff{P)=Tr[e-P"^"] 



N 

fc=i 



(14) 



Via Eq.(|^ one obtains the corresponding av- 
erage energy Ueff and the specific heat Ce//. 
Fig.[l] shows a plot of the average energy, com- 
paring the exact result with that from the effec- 
tive Hamiltonian. One observes that the agree- 
ment is better where T — > 0, i.e. in the low 
energy regime. A similar behavior is found for 
the specific heat, shown in Fig. [2]. In the limit 
P ^ CO the average energy tends to the ground 
state energy, U ^ huj/2 (Feynman-Kac formula) . 
One should keep in mind that H^f f has been con- 
structed for a specific value of the time parameter, 
T = 1 corresponding to the temperature T — I 
(we use ?i = fc_B = 1). The effective Hamiltonian, 
constructed in this way describes well thermody- 
namic observables in the range (3c < (3 (it works 
also for (3 > 10, not shown in the figure). How- 
ever, it breaks down for (3 < f3c, i-c. T > %. 
This is due to the small dimension = 20 of the 
matrix. Agreement in a larger /3-region, i.e. low- 
ering f3c can be obtained by increasing N. The 
statistical errors are of the size of the symbols in 
the Figs. Here we have not presented any results 
using the stochastic basis. This is discussed in 
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Harmonic oscillator 
m = l,h = l,uj = 0.6, Ax = 1, N = 20 
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Fig. [4] 
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Harmonic oscillator 



m = l, h = l, uj = 0.6, Ax = 1, N = 20 
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Fig. [5] 
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